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Abatract 

This paper presents a summary of some of the 
work performed by McDonnell Douglas Helicopter 
Company under NASA Langley-sponsored rotorcraft 
structural dynamics program known as DAMVIBS 
(Design Analysis Methods for VIBrationS). A set of 
guidelines which is applicable to dynamic modeling, 
analysis, testing, and correlation of both helicopter 
airframes and a large variety of structural finite ele- 
ment models is presented. Utilisation of these guide- 
lines and the key features of their applications to vi- 
bration modeling of helicopter airframes are discussed. 
Correlation studies with the test data, together with 
the development and applications of a set of efficient 
finite element model checkout procedures, are demon- 
strated on a large helicopter airframe finite element 
model. Finally, the lessons learned and the benefits 
resulting from this program are summarised. 

1. Background and Introduction 

Under a U. S. Government sponsored contract, Mc- 
Donnell Douglas Helicopter (MDHC) was one of the 
four helicopter manufacturers who participated in 
conducting a study of finite element analysis of he- 
licopter airframes for the enhancement of the tech- 
nology base related to the area of airframe structural 
dynamic analysis. This work which was sponsored by 
the NASA Langley Research Center Structures Direc- 
torate, consisted of planning, development, and the 
documentation of modeling techniques and computa- 
tional procedures. Additional tests were performed to 
obtain the necessary data to verify the finite element 
model and to develop procedures for correlation stud- 
ies. The modeling and shake test were performed on 
the McDonnell Douglas AH-04A (Apache) PV01 Heli- 
copter. Initially, a structural finite element model was 
prepared ***. 

« Mnnagar, Membsr Ali&. 

t Rsssareh Engineer, Member AHS. 

j Research Engineer, Member AHS. 


The guidelines were broken into three different cate- 
gories which consisted of those used for modeling of 
the stiffness characteristics, representation of the mass 
distribution, and consideration of factors affecting the 
accurate representation of vibrational characteristics 
of the model. Other resources required for performing 
efficient model checkout (2 \ and techniques used for 
different types of vibration analyses, were also devel- 
oped * 3,4 *. The portion of this effort which is related 
to the dynamic aspect of these studies is the subject 
of this paper. 

In the first section of this paper, the procedures used 
in converting the static model into a dynamic model 
are discussed. These include steps such as addition of 
non-structural components and extra grid points for 
the representation of the locations of the heavy mass 
items and inclusion of additional effects in shear pan- 
els to make them properly effective. For the accurate 
representation of the mass distribution, which is a key 
factor in a dynamic analysis, a procedure is developed 
which is discussed next ,6) . This procedure is used 
to systematically distribute the mass data over entire 
finite element model. In addition to these changes, a 
set of model quality assurance procedures is developed 
which helps the analyst to rapidly identify modeling 
errors prior to performing any analysis. These proce- 
dures are summarised. 

In the second part, the efforts spent on obtaining 
an experimental database describing the vibrational 
characteristics of the aircraft structure and the finite 
element model enhancement process are discussed * 
The procedure used for testing, the criteria used for 
mode identification, and a summary of test results 
is discussed. A summary of correlation studies per- 
formed to check the validity of a finite element model 
and the associated modeling assumptions will be dis- 
cussed. As a part of this enhancement process, the 
modeling problems that were identified as a result of 
correlation studies are presented. In addition, the gen- 
eral application of the test results is discussed. 

Finally, use of a model reduction procedure for the 
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purpose of increasing the computational efficiency for 
the dynamics optimisation applications is discussed. 
A preliminary study was performed to demonstrate 
application of such a procedure. The lessons learned 
from this study which include the development and 
application of the model reduction process and the use 
of the reduced model in a sensitivity analysis study 
are summarised. Finally the lessons learned from the 
overall DAMVIBS effort are discussed. 

11. Vibration Modeling 

Generally, there are two major tasks involved in the 
development of a dynamic model from a static model. 
First step involves the additional modeling resources 
necessary to make a static model useful for vibration 
analysis. Next is the generation and addition of the 
mass distribution to the static model. 

A static model which is primarily used for calcula- 
tions of internal loads and stresses under static loads 
generally does not contain a detailed representation of 
mass information. Prior to its use for dynamic anal- 
ysis, it is required that additional steps be taken to 
accommodate for the locations of heavy mass items 
such as engines, missiles, etc. In the case of Apache 
model, shown in Fig. 1, one of these steps involved the 
addition of secondary structures and non-structural 
parts and a set of additional grid locations required 
for proper representation of masses. These grid points 
were properly linked to the 



Fig. 1. AH-64A Dynamic MSC/NASTRAN Model. 

surrounding areas. The next step in this conversion 
process was to correctly represent the structure for the 
in service conditions. In the case of static models, gen* 
erally impending failure conditions are considered. In 
this case it is assumed that the skins are fully buckled 
and do not carry any tension/compression loads. This 
assumption needs to be modified for simulation of in 
service conditions. Thus, proper steps were taken in 
order to make sure that the skins were appropriately 
effective. 


Another task in generating a dynamic model involves 
the development of the mass model. This process gen- 
erally consists of following tasks; a) generation of a de- 
tailed weights record for the weight empty flight con- 
figuration, b) generation of a detailed listing of use- 
ful load weights, c.g.’s, and inertias, c) distribution of 
primary structure weight via material density parame- 
ter, d) manual distribution of large concentrated mass 
items (i.e., main rotor, transmission, engines, etc.), 
and e) automatic distribution of remaining mass data 
using a systematic and preferrably automated mass 
lumping process. 

Most of the information specified in the above list is 
generally available through the mass database and the 
mass density information of different geometric repre- 
sentation of various model components. For the pur- 
pose of distributing the the remaining mass data (i.e., 
item e), an automated mass lumping procedure was 
established. In the case of Apache model, this proce- 
dure was used to systematically distribute the mass 
data over the specified dynamic model grid locations 
and as a result a set of CONM2 cards were generated 
that were used in the dynamic model. The method- 
ology used for the distribution of the finite element 
model mass data is discussed in detail in Ref. 5. 

For the distribution of mass data, the lumping pro- 
gram requires two sets of information, namely, the ac- 
tual aircraft mass data records and the location of the 
selected grid points of the finite element model. The 
mass data is organised and stored in a format that is 
consistent with the M1L-STD1374A. Starting with the 
first mass item in the database, the program internally 
generates an imaginary volume, referred to as lump- 
ing volume, around the mass item. Next, by searching 
through the specified grid points, it identifies those 
grids which are confined within this volume and then 
assigns a different portion of the mass item to each of 
the selected grid points. Then, in case where there are 
still some remaining portions, the program increases 
the si se of the lumping volume and starts assigning 
portions of the remaining mass item to the new set 
of grids which may now be within the confinement of 
the increased volume. This process is repeated until 
the mass item is completely distributed to the sur- 
rounding grids. This process is repeated for each of 
the mass items. However, proper care is taken within 
the program for avoiding the relumping of the mass of 
those items which have already been accounted for. 

Figure 2 represents the general arrangement of the 
lumping process and illustrates the input and out- 
put streams of the automated mass lumping program. 
Although this program can accommodate structural 
mass input, the primary intent in the development 
of the program was to distribute nonstructural mass 
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items («.«• fasteners, wire bundles, etc.) to model 
grid point locations. This mass lumping process min- 
imises the human error through automatic generation 
of a set of C0NM2 cards which are readily used as 



Fig. 2. Block Diagram of Mass Lumping Program. 


a part of the NASTRAN bulk data deck. The struc- 
tural mass (e.g. skins and stringers) is calculated via 
material density cards internally within NASTRAN 
to form the total mass matrix. 

IIL Model Yerificfttioa .Proew 

Model verification process is one of many tasks re- 
quiring a substantial amount of time. Thus, it is very 
desirable to have tools that can facilitate and auto- 
mate this process as much as possible. A series of 
computational procedures were developed in the form 
of DMAPs. These DMAPs were used to identify both 
modeling errors and also provide certain information 
regarding identification of modes, energy associated 
with each mode and calculation of modal participation 
factors. The latter two DMAPs were developed prior 
to start of DAMVIBS program at the MDHC. How- 
ever, due to their high level of effectiveness and also 
for the sake of completeness they are included here 
and are briefly discussed in the following sections. 

As a result of applications of these modeling check- 
out DMAPs to a variety of structures, a common set 
of sources of errors has been identified. Those errors 
which occur during the assembly of the G-set stiffness 
matrix are due to; (a) the improper specification of 
the coordinate couplings where two components which 
are located in different coordinate systems are cou- 
pled together improperly, (f>) the use of a short beam 
element next to a long beam element, (c) use of a 
large or improper aspect ratio for plate elements, or 
(d) the use of CELAS elements between noncoinci- 
dent points. Each of these modeling practices will 
result in the ill-conditioning of the G-set stiffness ma- 
trix. In the second and third levels of model forma- _ 


tion, where the N-set and F-set stiffness matrices are 
being assembled t improper applications of the Multi- 
point Constraint equations (MPCs) and Single Point 
Constraints (SPCs) could result in further modeling 
errors. Incorrect specifications of these types of con- 
straints will result in incorrect representation of the 
linear relationships among different DOFs and over- 
constrained boundary conditions, respectively. 

C holesky Decomposition DMAP 

The purpose of this check is to identify singularities or 
mechanisms as well as near singularities or soft spots 
in the model at the F-set level. For this purpose, a 
DMAP was developed which employs the DCOMP 
functional module to perform a triangular decomposi- 
tion of the F-set stiffness matrix in the static solution 
sequence [2) . After the decomposition, a diagonal ma- 
trix results whose elements provide information about 
the conditioning of the stiffness matrix. A max factor 
diagonal value (i.e., ratio of the largest diagonal term 
to the lowest diagonal term) of over 10 6 indicates a soft 
spot (i.e., near singularity) while a value over 10 10 is 
indicative of a mechanism (i.e., singularity). In ad- 
dition to the max factor ratio, another parameter, 6, 
is calculated. Aside from round-off errors, any large 
nonzero value of € is indicative of a singularity in the 
stiffness matrix. 

Multi-Level Strain Energy DMA P 

In contrast to the Choleski Decomposition Check, 
Multi-Level Strain Energy Check provides informa- 
tion about modeling problems at all three levels (i.e., 
G-set, N-set, and F-set) of model formation in NAS- 
TRAN. Through a DMAP, the stiffness matrix, and 
the rigid body modes obtained from the grid point 
geometry of the structure, are used to calculate the 
strain energy of the structure at each of the three 
levels of model formation. The same information is 
also used to calculate another parameter referred to as 
the nodal strain to provide further information about 
modeling problems at each of the levels |21 . As a 
result, any problems which might be caused by the 
improper grounding of the structure (i.e., G-set), in- 
correct application of MPC equations (i.e., N-set), or 
incorrect definition of the SPC constraints (i.e., F-set) 
will be identified through examination of these nodal 
strain energy at each level. To demonstrate the ap- 
plicability of this powerful tool, it was applied to the 
AH-1G NASTRAN model Examination of the N- 
set strain energy matrix, shown in Table 1, revealed 
some problems associated with the rotational degrees 
of freedoms (i.e., the third through sixth diagonal el- 
ements of the matrix are much larger than 10) which 
are indicative of incorrect application of MPC equa- 
tions. To further identify the degree of freedom as- 
sociated with the MPC equation, another matrix, re- 
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ferred to u N-iet nodal strain matrix was examined. 
The last 
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Table l. N-set Strain Energy Matrix. 


three columns of this matrix, shown in Table 2, which 
were representative of the rotational degrees of free- 
dom revealed the location of the troublesome grid and 
the associated degrees of freedom (i*e., grid points 
7505 and 15218). Once the associated MPC equation 
was corrected the DMAP was again applied and sim- 
ilar problem associated with the MPC equation used 
for the engine support location was identified. This 
problem was then corrected. 
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Table 2. N-set Nodal Strain Matrix. 

Kinetic Energy DMAP 

One of the key recent improvements in the dynamic 
analysis of large finite element models is the devel- 
opment of efficient and accurate eigenvalue extraction 
techniques. Introduction of such techniques has re- 
sulted in their applications to more complex systems 
and consequently much larger FEMs. As a result, for 
the prediction of structural modal characteristics, the 
analyst is often confronted with the identification and 
post processing of large number of high density modal 
vectors. For very large and detailed models, this pro- 
cess is rather tedious and time consuming. Conse- 
quently a DMAP was developed which is referred to 
as Kinetic Energy Check DMAP ,3,4) . Application of 
this DMAP results in the caculation of the modal ki- 
netic energy associated with each mode and rapidly 
identifies the candidate structural modes. Addition- 
ally, such results provide the necessary information to 
identify the local modes and any possible modeling er- 
rors. Since the development of this DMAP, it has been 
applied to a large number of different structures and 
has resulted in a substantial saving of manhours and 
the elimination of costly and time consuming plotting 
process of the local/trivial modes. 


This technique was initially applied during dynamic 
analysis of the Hughes Advanced Attack Helicopter 
(A AH). At the time of application, the model had 9792 
degrees of freedom. To identify the modal characteris- 
tics of the model, an eigen analysis was performed up 
to 25 Hs. Within this frequency range, eighty seven 
modes were calculted where each mode contained an 
eigenvector containing 10000 values. Through appli- 
cation of the kinetic energy DMAP, a filtered modal 
kinetic energy matrix was obtained which contained 
only a set of larger values for each mode ,4> . A quick 
examination of the filtered modal kinetic energy ma- 
trix revealed candidates for the local or actual struc- 
tural modes. A structural mode is generally identified 
by a column containing many large numbers, since the 
energy is distributed and large portion of structure 
moves together. In the case of local modes, only a few 
elements of the matrix have very large values, indicat- 
ing only a local area of structure is moving. Figure 3 
shows a sample of a structural mode identified through 
this procedure. 



Fig. 3. AAH First Torsion Mode. 


Mod a l Participation DMA P 

An important information during studies of the dy- 
namic structural response of linear systems is the con- 
tribution of different structural modes to the total re- 
sponse of a selected point or a finite region of the 
structure. Such information generally provides valu- 
able insight into the vibrational characteristics of the 
system and points out areas where structural modifi- 
cations could result in improved response. In certin 
applications, such data provides the necessary infor- 
mation for the identification and selection of an opti- 
mum location for vibration control devices. For this 
purpose, a DMAP was developed which calculates the 
modal participation factors associated with response 
of a selected location of the structure For the ap- 
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plication of this DMAF, frequency response (SOL 30) 
is used. Application of this technique to a vibration 
reduction study is discussed in detail in Ref. 5. 


IV. Experimental Efforts 


Test Objective 

The main objective of the shake test was to obtain an 
experimental data base describing the vibration re- 
sponse of the test vehicle. Thu data base consuted 
mainly of transfer functions of response acceleration 
to input force versus input frequency for various mea- 
surement locations throughout the fuselage structure. 
Although the test was not a modal survey test, the 
major airframe modes were extracted from the test 
data. These modes provided insight into the basic vi- 
bration characteristics of the fuselage and aided in the 
correlation with NASTRAN model. 

Test Vehicle 

The test vehicle used was the U.S. Army’s AH-64A 
Advanced Attack Helicopter. The testing was per- 
formed with the vehicle in the primary mission con- 
figuration ,<M . The main rotor hub and blades were 
replaced with a rigid steel fixture which was attached 
to the static mast in a manner similar to the actual 
hub. The dummy structure was equivalent to 100% 
of the weight of all of the actual items except for the 
rotor blades. The blade weight was reduced to 60% 
of the actual weight to more closely represent the dy- 
namic equivalent of the rotating system. 

Test Setup 

The test vehicle was suspended from a gantry. Air sus- 
pension springs located at the top of the gantry pro- 
vided vibration isolation between the test vehicle and 
the gantry structure. This method minimised feed- 
back from the gantry modes of vibration into the test 
aircraft. Link chain was used to attach the air spring 
system to the dummy main rotor hub. This technique 
removed the tendency for the aircraft to yaw. Sinu- 
soidal vibration excitation was applied to the test ve- 
hicle at the main rotor hub and at the tail rotor. Three 
forces (longitudinal, lateral and vertical) were applied 
to the main rotor hub, as well as two moments (pitch 
and roll). These forces and moments were applied 
at the vertical and station positions equivalent to the 
blade plane and the center of rotation. Three forces 
were applied at the tail rotor location: longitudinal, 
lateral and vertical. The tail rotor excitation point 
was at the center of rotation and at a lateral position 
which is half way between the two rotor planes. The 
flight tail rotor equipment was on the aircraft during 
the entire shake test. 


Test Procedure 

During excitation, data was taken over a frequency 
range of 3.75 to 50 Hs except for the tail rotor in- 
put cases which extended this range to 200 Hs. The 
frequency was stepped up in increments of 0.5% and 
the input dwelled for a few seconds at each frequency 
step to allow the response to come to quasi-steady 
state. Data was taken for all accelerometers at each 
frequency step. Each dwell and measurement step 
took about 5 seconds to complete, so an entire test 
(3.75 to 50 Hs) could be run in less than 45 minutes. 

Instrumentation 

A total of 102 accelerometers were mounted on the 
aircraft at 39 locations and four were mounted on the 
gantry and suspension system. Each accelerometer 
was tested in a mechanical system, comprised of a 
load cell and a known free mass. The mass was driven 
in the test frequency range. The resulting transfer 
function from each unit was used by the computer to 
convert voltage to acceleration and to adjust the phase 
at each frequency step. 

The accelerometer locations were chosen based on 
three criteria: (1) to identify the response at key lo- 
cations such as pilot seats, engines, hubs, etc., (2) to 
describe fundamental modes of the airframe; and (3) 
to identify important local modes such as those of the 
wings. 

The accelerometers were mounted on one inch cubic 
fiberglass blocks which were epoxied to the airframe. 
The blocks were oriented such that their faces were 
aligned with three orthogonal axes of the ship. Each 
block had one to three accelerometers attached to it. 

The use of a large number of fixed accelerometers 
proved to be much more effective than a few roving 
accelerometers, as is usually done in this type of test- 
ing. The chief reasons for this are as follows. The 
accuracy of accelerometer positioning was greatly in- 
creased. Errors in mounting and polarity were min- 
imised because the accelerometers remained fixed in 
position for the entire test. Moreover, the time re- 
quired for testing was greatly reduced because only 
one frequency sweep was necessary for each test case. 

Data Acq uisition and Control 

The system used for data acquisition and control 
was 100% computer controlled and exclusively used 
Hewlett Packard instrumentation. At the heart of the 
system was the HP9836 computer/system controller. 

At each frequency step, all data channels were scanned 
and 16 readings were taken over one cycle of vibration. 
This data was then sent from the HP6942A scanner 
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to the HP9836 computer where magnitude and phase 
data was extracted by means of a Discrete Fourier 
Transform. By this method the first harmonic was 
extracted while the higher harmonics, up to the 8th, 
were rejected. 



Fig. 4. Fuselage First Vertical Bending, 5.45 Hs, 
Test Results 

The essential results of the shake test were the re- 
sponse transfer functions. Additional results include 
natural frequencies and mode shapes which were 
gleaned from the transfer functions. Although not es- 
sential, this data aided in the correlation of the NAS- 
TRAN model and provided insight into the vibration 
characteristics of the aircraft. A sample mode shape 
from the test data is presented in Fig. 4. Natural 
frequenies were estimated from the data by means of 
an indicator function When the value of the first 
function approaches unity, a natural frequency is in- 
dicated. Sample indicator function results are shown 
in Fig. 5. 



Fig. 5. Sample Indicator Function Results. 

v. Correlation and Mo deLJEnhaiKcmcRt 

An excellent correspondence between test and NAS- 
TRAN results was obtained for the frequency re- 
sponses up to 13 Hs. This includes the first six modes 
found in the shake test. Above 13 Hs, the correla- 


tion tends to deteriorate gradually as the frequency of 
the response increases. A comparison of natural fre- 
quencies found by test and calculated by NASTRAN 
is given in Table 3. Typical response correlation plots 


are ihown in Fig- 6. 
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TAILBOCH TORSION 

4.42 

4.28 

-3.2 

1st VERTICAL BENDING 

5.45 

5.39 

-1.1 

1st LATERAL BENDING 

9.39 

0.79 

-6.4 

LONG. BENDING OF VERTICAL 

10.07 

9.78 

-2.9 

IVR MAST BENDING, ICNG. 

12.11 

12.26 

1.2 

m/r mast BENDING, LAX. 

12.48 

12.74 

2.1 

SYMMETRIC WING BENDING 

13.38 

15.15 

13.2 

Atrn-svm. WING BENDING 

13.85 

15.37 

11.0 

STABILATOR YAW 

20.44 

17.78 

-13.0 


Table 3. Correlated Modes 


The problem areas in the correlation were the en- 
gines, stabilator, and the wings. However, much in- 
provement was obtained in the vertical response of the 
wings based on the preliminary correlation results. In 
addition to improving the existing finite model, much 
knowledge gained in the correlation effort will be ap- 
plicable to future analyses. Shake test results can of- 
ten be used to refine an existing finite element model. 
This is accomplished by comparisons of experimental 
and analytical results. This comparison initially in- 
dicates whether or not there are errors in the model 
and assists in locating probable areas in which these 
errors may reside. If an error has been caused by an 
incorrect assumption, then the correlation study can 
provide a guide as to how the problem might be cor- 
rected. Finally, the test comparison provides a stan- 
dard by which the degree of enhancement can be mea- 
sured. The following discussion describes an example 
of how shake test results were utilised to enhance the 
NASTRAN finite element model of the AH-64A air- 
frame. 

Initial comparison of the frequency response curves 
obtained from the test with those predicted by NAS- 
TRAN indicated that, for some cases, there were dis- 
crepancies between the two sets of results. The most 
noticeable problems appeared in the responses of the 
wings, the engines and the stabliator. 

Examination of the deflected shapes obtained from 
test and analysis helped to further localise and iden- 
tify errors in the model. A specific example of such er- 
ror is a mode at approximately 8.5 Hs which is shown 
in Fig. 7. This mode is predominantly a wing mode 
although it is significantly lower in frequency than the 
wing modes found in the test results. In addition, one 
wing is undergoing pure vertical bending, while the 
other is experiencing nearly pure wing torsion. Such 
asymmetric behavior was not found in the test results. 
Therefore, it was concluded that there were errors in 
the wing portion of the finite element model. 
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After clou .crutiny of the wing portion of the model, 
problems were identified and corrected. The problems 
were mainly asymmetries in the lumped masses and 
element products of inertia, overly stiff modeling of 
the pylons, and the omission of the trailing edge. The 
wing trailing edge was ignored in the static model be- 
cause it was considered secondary structure. 



(a) Longitudinal Raaponaa at Main 
Rotor Hub 



(b) Vurtical Raaponaa at Aircraft 

Noaa 




(d) vartical Raaponaa at Engina 



Fif. 7. Wing Mode at 8.5 Hs. 
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With the improvements made to the model, the ana- 
lyst can use the entire inventory of test results such as 
frequency response plots, deflected shapes, and modal 
frequencies to examine the effect of the improvements. 
An excellent example of this is the improvement of the 
vertical response of the wing as a result of the changes 
implemented in the model. See Pig. 8. 



tIMNMQ KIM. KPMU 

Fig. 8. Improved Correlation of 

Vertical Response at Wing Tip. 

VI. General Applicatio n of T eat Results 

A comprehensive shake test, such a a the one per- 
formed under this program, not only helps to enhance 
an existing finite element model of the test specimen, 
but also provides valuable guides for enhancing other 
models and future modeling efforts. Test results may 
be used to verify modeling practices used in the past, 
to identify questionable assumptions and practices, 
and to provide a reasonable level of confidence in the 
analysis. 

It is usually necessary to make assumptions in prepar- 
ing a finite element model. The validity of these 
assumptions must be verified by testing. In this 
particular case, the test results confirmed many of 
the assumptions used in modeling the AH-64A heli- 
copter. For example, the assumption that stringers 
and longerons carry axial loads while skins and webs 
carry shear loads was verified. The effectiveness of 
using a static model for dynamic analysis with dy- 
namic reduction was also confirmed. Good correlation 
associated with main rotor mast responses confirmed 
our modeling technique regarding the main rotor sup- 
port structure and the overall representation of the 
airframe. 

The correlation of the test results revealed some as- 
sumptions and practices which were not quite correct. 
The following are some examples. Although using a 


static model for dynamic analysis was generally suc- 
cessful, some of the assumptions made for statics were 
not good for a dynamic model. Two examples were 
ignoring the trailing edge structure of the wing and 
the frictional load carrying capability of bolts. Some 
assumptions made for statics were previously recog- 
nised and corrected to meet dynamics requirements. 
An example of this was the shear panel effectivity. 
Regarding the mass distribution, some of the assump- 
tions and methods were not strictly correct. For ex- 
ample, the automated mass lumping routine does con- 
sider connectivity when lumping the masses, nor does 
it rigorously account for rotational mass moments of 
inertia. 

VIL Vibration Reduction Studies 

Generally, use of model reduction procedures for the 
purpose of increasing the computational efficiency is 
of great interest in a dynamic analysis. This is espe- 
cially true in dynamic optimisation applications. Con- 
sequently, a preliminary study was done to explore the 
applicability of a reduction procedure to a vibration 
reduction study. A model reduction procedure was 
initially developed and was used in vibration study. 


A computational procedure for the reduction of large 
airframe finite element models was developed ***. This 
procedure, which has been implemented as a set of 
DMAPs, results in a significantly reduced model while 
retaining the essential dynamic characteristics of the 
full-sised model. This procedure was applied to a 
slightly modified version of the airframe dynamic fi- 
nite element model of AH-64A Attack Helicopter. A 
mathematicaly reduced model with significantly less 
DOFs was obtained. This reduced model, is shown in 



This model is an accurate representation of the global 
behavior of the full model. The resulting reduced 
model resulted in a substantial reduction in the com- 
putation time. This reduced model which has 83 grid 
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points and a total of 498 DOF is only a mathematical 
representation of the reduced (stick) model (i.e., the 
elastic properties are defined in terms of stiffness ma- 
trix rather than physical elements such as bar, beam, 
etc.)- However, further steps could be taken to carry 
this mathematical reduced model one step further and 
develop an equivalent “physical” reduced model. This 
process was not pursued. Instead, the “mathemati- 
cal” reduced model was adapted for sensitivity analy- 
sis study which will be discussed in the following sec- 
tion. Following the reduction process, an eigen anal- 
ysis was performed to evaluae the accuracy of the dy- 
namic reduced model. Figure 10, represents the first 
vertical bending mode. 




Fig. 10. Stick Model First Vertical Bending. 

To further quantify the degree of correlation between 
the two models, an in-house computer program which 
provided a systematic way of checking the degree of 
correlation between frequencies, mode shapes or a 
combination of the two, was used (8) . Table 4 shows 
the resulting correlation matrix between the two mod- 
els. 

FREQUENCY (HZ) MODE SHAPE 

MODe full stic k corr elation* 

TAILBOOM TORSiON * 45 

1ST VERT BENDING «00 

1ST LATERAL BENDING 10 70 

SYM ENGINE YAW AND PITCH 11 44 

VERT TAIL LONG BENDING II 97 

MAST LONG. BENDING IXH 

ANTISYM ENGINE YAW 14 16 

8TABILATOR YAW 


Stt 

on 

6 II 

093 

9 76 

064 

II 67 

090 

11.31 

091 

14 33 

0 67 

16 43 

0 70 

I960 

OBI 


* 1.00 s PERFECT CORRELATION 
0 00 - NO CORRELATION 

Table 4. Modal Correlation Matrix. 

The correlation factors indicated in Table 4 are basi- 
cally indicative of the degree of correlation between a 
pair of global modes. 

Vibration Reduction study 


In addition to the development and application of the 
model reduction procedure, another study was made 
to examine the applicability of the reduced (stick) 
model to a vibration reduction study. For this pur- 
pose, the reduced model was subjected to different 
four-per-rev hub excitations and the modal frequency 
responses (SOL 30) of different locations of the struc- 
ture, together with the contribution of each mode to 


the total response of each location, were calculated 
l 7 *. Subsequent to the identification of the dominant 
modes, a design sensitivity analysis study (SOL 53) 
was performed to identify the pertinent model param- 
eters (e.g., cross sectional area, area moment of iner- 
tia, etc.) which have the most effect on the vibrational 
response at the selected locations. Once these param- 
eters were identified, certain incremental changes were 
made to each of them individually or a combination 
of these parameters to reduce the vibration at the se- 
lected locations. 


In conjunction with the vibration reduction study, the 
contribution of each mode to the total response of the 
selected portion of the structure was determined at the 
location of interest. This location was the tip of the 
vertical tail of the AH-64A reduced model. For this 
purpose, an in-house DMAP (7) program was used 
in conjunction with the MSC/NASTRAN modal fre- 
quency response solution sequence (SOL 30). Figure 
11 shows the four modes with the largest contribu- 
tions to the total response of the selected point, when 
the aircraft is subjected to a four-per-rev vertical hub 
excitation. 



STICK MODEL 

MO0C-I7. 10 603 HZ 
HO0C-2O. 21900 HZ 
MOOf 0. 6 153 HI 


- MODS- 7. 5 623 HZ 

MOOC-15. 17 266 HZ 

TOTAL MSFONSt 


Fig. 11. Structural Modal Response. 


As can be seen, mode 17 contributes the most to 
the total response. In addition, examination of other 
modes and the results obtained from the DMAP in- 
dicated that the fifteenth mode (i.e., 17.268 Hs) to 
be the next highest contributor to the response of the 
point of interest. Therefore, these two modes (i.e., 
modes 15 and 17) were selected to be used in a follow- 
ing design sensitivity analysis. This selection process 
was repeated for two other types of hub excitations 
(i.e., longitudinal and lateral excitations) and a simi- 
lar pattern was obtained. 

To identify those parameters which have the most ef- 
fect on the response of selected normal modes in the 
area of interest of the structure, a sensitivity analysis 
was performed. For this study, the tailboom which 
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is in the close proximity of the vertical tail was se- 
lected. The cross-sectional area, the two area mo- 
ments of inerti& y and the torsional constant parame- 
ters were selected as the design variables. The two 
dominant stick model modes (i.e., modes 15 and 17) 
which represented mainly the tailboom vertical and 
lateral bendings, were selected as the “constraint” pa- 
rameters. The design sensitivity analysis was per- 
formed for all twelve tailboom frame segements and 
a set of sensitivity coefficients obtained. These are 
shown in Fig. 11 and 12. 
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Fig. 11. Sensitivity Coefficients, Mode 15 
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Fig. 12. Sensitivity Coefficients, Mode 17 

Examination of these figures revealed that the two 
area moments of inertia parameters have the most ef- 
fect on the frequency placement of the two selected 
modes. Consequently, these parameters were altered 
which shifted the frequencies of these two modes away 
from the four-per-rev excitation frequency < 91 . This 
process resulted in a reduction of the vibration level 
at the tip of vertical tail. Similar changes were also 
made to the same locations of the full model and the 
response of the vertical tail was calculated. 


VIII. Concluding Remarks 


A summary of the MDHC’s effort related to the devel- 
opment of DAMVIBS dynamic modeling guidelines, 
testing and the correlation studies has been presented. 
As a result of this work, a set of comprehensive mod- 
eling procedure and check out tools have been gener- 
ated which have been used in modeling of large air- 
frame structures. These procedures which been used 
in the ongoing enhancement process of the exiting 
MDHC product models, have also been applied to the 
newly developed product models. Throughout this 
program, a substantial amount of experienced have 
been gained in creating other reliable airframe finite 
element models (10) which are used during the design 
process. Through such improved modeling capability, 
a better estimation for man-hour requirements and 
scheduling process is now possible. 

In relation to performing shake tests and efficient use 
of this data, our capability has been substantially en- 
hanced. This has been accomplished through use of 
experience gained during this program in the preper- 
ation of better test plans, performing the the shake 
tests where issues such as support of test structure, 
applications of loads or load levels are of great im- 
portance. Reduction and application of test data to 
the model resulted in improved corrolation method- 
ology which also brought into focus the shortcomings 
associated with different correlation methods and the 
inaccuracy associated with higher analytical modes. 
Further examination of the test data pointed out the 
degree of structural nonlinearity and provided insight 
into damping characteristics of the aircraft structure. 

The exchange of information between MDHC and 
other three companies provided a means to compare 
methodologies used by each group and the associated 
results. Such an exchange of information resulted in 
improving certain areas of our models. One of these 
areas was inclusion of secondary structural compo- 
nents which have resulted in improving the correla- 
tion results. Finall, these studies have pointed out 
new challenges in terms of limitations associated with 
FEA which need to be further investigated. 

A summary of the lessons learned during these studies 
are listed as follow: 

• The static finite element model may be used for 
dynamic analysis after making the proper mod- 
ifications. 

• The proper modeling of the mass distribu- 
tion and representation of secondary and uon- 
structural components is essential to accurate 
vibration modeling. Automation of this proce- 
dure greatly reduces modeling time. 
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• Model checkout and verification facilitates ob- 
taining accurate results. Automated DMAP 
procedures enhance the analysts capabilities in 
this area. 

• The performance of the finite element model can 
be significantly improved by the application of 
the shake test results. 

• Utilization of stationary accelerometers reduces 
the test time. 

• Prior to performing the test, magnitude and 
types of excitation loads should be studied using 
FEM. 

• Excellent correlation with AH-64 shake test re- 
sults was obtained at most measurement loca- 
tions up to 13 Hz. 

• Knowledge gained in such a correlation study is 
also applicable to future analysis efforts. 

• The vibration reduction obtained as a result 
of a preliminary dynamic optimization study 
demonstrated the applicability and benefits of 
the model reduction technique. 

• Confidence in finite element analysis as a valid 
tool for predicting helicopter fuselage vibrations 
has been greatly increased. 
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